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ABSTRACT 

We predict how protoplanetary disks around low-mass young stars would 
appear in molecular lines observed with the ALMA interferometer. Our goal 
is to identify those molecules and transitions that can be used to probe and 
distinguish between chemical and physical disk structure and to define neces- 
sary requirements for ALMA observations. Disk models with and without verti- 
cal temperature gradient as well as with uniform abundances and those from a 
chemical network are considered. As an example, we show the channel maps of 
HCO"'"(4-3) synthesized with a non-LTE line radiative transfer code and used as 
an input to the GILDAS ALMA simulator to produce noise-added realistic im- 
ages. The channel maps reveal complex asymmetric patterns even for the model 
with uniform abundances and no vertical thermal gradient. We find that a spa- 
tial resolution of 0.2 — 0.5" and 0.5-10 hours of integration time will be needed 
to disentangle large-scale temperature gradients and the chemical stratification 
in disks in lines of abundant molecules. 

Subject headings: accretion, accretion disks — astrochemistry — circumstellar 
matter — molecular processes — planetary systems: protoplanetary disks — 
instrumentation: interferometers 



1. Introduction 

The cradles of planet formation - protoplanetary disks - are ubiquitous around young 
stars, at least of low and intermediate masses (< 8 Mq). It is imperative to study disk struc- 
ture and evolution in detail if one aims to understand initial conditions for planet formation. 
These objects are best studied at (sub-) millimeter wavelengths in thermal dust continuum 
and various molecular lines tracing their physical and chemical structure. However, being 
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small (~ 200 — 1000 AU) and of relatively low mass (~ O.OIM©), such disks are hard to 
study observationally as high spatial resolution and sensitivity are required. To date only 
several nearby disks around T Tauri and Herbig A e/Be stars have been spatially resolved 
in a handful of molecular lines (IBergin et al.l 120071 ) . These observational data coupled to 
advanced modeling allowed us to constrain basic disk parameters such as size, mass, temper 



ature, density distribution, kinematics, and ionizat ion structure (e.g.. ISemenov et al.l 12005 
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tical temperature gradient that seems to be present iri some disks (e.g., DM Tau, AB Aur) 



while absent in others (e.g., LkCa 15), see iPietu et al.l (l2005l . 120071 ) 



The situation will change dramatically when the Atacama Large Millimeter Array 
(ALMA), equipped with 50 12-m antennas, will come into operation in 2012. ALMA will 
be capable of imaging protoplanetary disks at spatial resolutions up to ~ 0.04 — 0.005" in 
a frequency range between 100 and 950 GHz. This will allow direct detection and char- 
acterization of disk instabilities, resulting in clumpy structures (vorticit ies, "spiral arms" , 
etc.) as well as inner gaps and holes induced by forming giant pla nets (IWolf et al.ll2002al : 
Wolf fc Klalii]l2002bl : Iwolf fc D'Angeloll200,4 iNaravanan et aPboOfil ). The large-scale chem- 
ical structure of pro toplanetary disks will beco r ne accessible in rotat ional lines of many 
ab undant molecules (lAsensio Ramos et al.l 120071 ; iKrumholz et al.l l2007l ) . As demonstrated 
by iPavlyuchenkov et al.l (120071 ) . excitation of molecular lines in disks with strong gradients 
of physical conditions and chemical structure can be a complicated process and may be hard 
to interpret. 

The main aim of this work is to unravel the potential of ALMA to distinguish between 
various temperature and chemical effects in protoplanetary disks, with an emphasis on dis- 
cerning spatial resolution and integration time needed for this goal. Such observations will 
allow us to fully explore the validity of state-of-the-art disk chemical models on a global 
scale for the first time, and a dvance our unders t andin g of the disk physics enormously. In 
contrast to previous studies of iNarayanan et al.l (120061 ) and lKrumholz et al.l (120071 ). we take 
into account both the non-LTE excitation and realistic chemical structure. 



2. Disk models 



We adopt a steady-state flared disk model with a radius of 800 AU, an accretion 
rate M = 10~^ M^r^ y r~^, a viscosity parameter a = 0.01, and a mass of M = 0.07 



(iD'Alessio et al.lll999l ). The disk has a microturbulent velocity of 0.1 kms ^ and a Keple- 

This disk model rese mbles the disk arou nd the young 
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rian rotation profile, V{r) 
T Tauri star DM Tau located at the distance of ~ 140 pc (iDutrey et al.l 120071 ). Its thermal 
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and density structure is shown in Fig. [T] (left and middle panels). To mimic the case when 
the vertical temperature gradient is weak or absent, we also consider the same disk model, 
but with a vertical temperature distribution fixed to the value at one pressure scale height. 

Using both disk models, 5 Myr o f disk chemical evolution is modeled with a gas-grain 
chemical model from Vasvunin et al. (boOTh. The gas- phase reaction data are adopted from 
the RATE 06 database (jWoodall et al.l 120071 ). Ionization and dissociation of the gas-phase 
species as well as desorption of surface species due to cosmic rays, stellar X-rays, stellar 
and interstellar UV radiation are taken into account. The disk is illuminated by UV radia- 
tion from the central star with ar i intens i ty G = 410 Go at r = 100 AU and by interstellar 
UV radiation with intensity Gq (jPraind 1 19781 : Ivan Dishoeckl 119881). The X -ray ionization 



rate in a given disk region is computed acco rding to iGlassgold et al.l (119971) . The adopted 



cosmic-ray ionization rate is 1.3x10 s (ISpitzer &: Tomaskol Il968l ). Molecules are as- 



sumed to stick to grains with 100% efficiency. We also consider dissociative recombination 
of ions on charged dust grains and grain re-charging, with a dust-to-gas mass ratio of 1%. 
To model surface chemistry on th e 0. 1/xm olivine grains, a set o f surface reactions from 
Hasegawa. Herbst. fc Leung) (1l992l ) and iHasegawa fc HerbstI (1l993l ) is used. 



In what follows, we will concentrate on HCO^ - a representative abundant molecule that 
traces the ionizatio n fraction and posse sses strong rotational transitions, and which is readily 
observed in disks (IDutrey et al.ll2007l ). The modeled distribution of the HCO"*" absolute 
abundances (molecules cm~^) in the DM Tau disk at 5 Myr is shown in Fig. [T] (right panel). 
A layered chemical structure with the maximal HCO"'" concentration reached at intermediate 
disk heights is c learly visible. Many ob s erved species have a similar stratification, e.g., CO, 
CS, HCN, etc. kikawa fc HerbstI 1 19991 : fvasvunin et al.ll2007h . 



Strong turbulent mixing and/or global advection fio ws can smooth these abundance gra- 
dient s, leading to more uniformly distributed abundances (IWillacy et al. I bood : ISemenov et al.l 
20061 ). Therefore, we also consider a model with an uniform HCO^ abundance of 10~^ rel- 
ative to H2. Thus, our study is based on three disk models: (1) the model with chemical 
stratification and vertical temperature gradient, (2) the same model but with uniform abun- 
dances, and (3) the model with uniform abundances and no vertical temperature gradient. 



3. Synthetic channel maps 



Using these three DM Tau-l i ke dis k models, and the 2D non-LTE line radiative trans- 
fer code of iPavlyuchenkov et al.l (120071 ) with thermal dust continuum included, we synthe- 
size "ideal" (beam-unconvolved) channel maps. The coUisional rates are taken from the 
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Leiden Atomic and Molecular Databasjl 



(ISchoier et al.ll2005l ). while dust opacities for sili- 



cate grains are from iDraine fc Led (Il984j ). At the 60° inclination the integrated HC0"'"(4- 



3) line wi dth is about 3 kms \ wit h a peak intensity reached at about 0.8 kms ^ (see 



Fig. 10 in iPavlyuchenkov et al.l 120071 ). The corresponding 0.1 kms ^ wide velocity channel 



at V = 0.77 kms ^ of the continuum-subtracted HCO"*"(4-3) map is shown in Fig. [2l 

This close to edge-on orientation is particularly favorable as the layered and uniform 
disk chemical structures can be clearly distinguished in the channel maps, if sufficiently small 
channel widths of 0.1 — 0.2 kms~^ are used. All synthetic channel maps reveal a complex 
pattern and are asymmetric. The model with both non-zero temperature and abundance 
gradients has a remarkable "omega" shape, where the cold midplane with low HCO"'" concen- 
tration appears as two intensity "holes" (Fig. [2], left panel) . A steep temperature gradient in 
radial direction toward the inner disk region is clearly visible, while the vertical temperature 
gradient cannot be fully traced with the HCO^ lines due to strong chemical stratification of 
this species. The HCO"'" (4-3) emission is thermalized and optically thin in this case. 

The two models with uniform abundances have relatively high HCO+ column densities 
so that self-absorption becomes important. In the model with fixed vertical temperature this 
effect leads to a layer with low intensity in the upper part of the map (Fig. [21 right panel). 
However, the overall intensity of the HCO"*" (4-3) emission is higher compared to the model 
with layered abundances. The radial temperature gradient is also prominent in this case. 

The model without chemical stratification but with vertical temperature gradient also 
reveals the self- absorption layer in the upper part of the channel map. The excitation 
temperature of the HCO"*" (4-3) transition changes strongly with the disk height and has the 
lowest value in the midplane. This zone of low intensity appears as a fake "spiral arm" in 
the lower part of the map (marked as "midplane" in Fig. [21 middle panel). Note that such a 
pattern somewhat resembles the "omega" structure for the model with chemical stratification 
of HCO"*". While we show only one representative channel, one should bear in mind that 
analysis of the interferometric data is more reliable when all channels and dust continuum 
emission are taken into account. 



4. ALMA observations 

Despite the fact that all the models considered can be distinguished in the synthetic 
channel maps, it is important to study whether ALMA is able to disentangle the effects 
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of thermal and chemical structure in disks. To simulate the observations with ALMA, we 
us e the synth e tic H CQ"'"(4-3) channel maps as an input to the GILDAS simulation software 



of iPety et al.l (1200 ll ) (version from Oct ober 2007). A positiori close to DM Tau is assumed 



4h 33m 48s, ^2000 = 18° 10' 09" JPucourant et al.lboosh . We focus on the 0.77 km 



s-i 



0^2000 

channel (with 120 kHz bandwidth or 0.1 kms~^ velocity width). Array configurations with 
64 antennas built in GILDAS are considered. Note that in reality ALMA will consist of 50 
antennas with different configurations, so our analysis will tend to underestimate required 
observing time by a factor of 1.5 but the simulated spectra will not be significantly affected. 



Typical weather conditions at the Chajnantor plateau are assumed (see iPety et al. 



2OOII ). with the following main types of errors leading to n oise: 1) receiver temperature 



is 80 K, system temperature is 230 K {at u = 357 GHz, see lGuilloteaijl2002l ). 2) random 
pointing errors during the observation are 0.6", 3) relative amplitude errors are 3% with a 
6% hour~^ drift, 4) residual phase noise after calibration is 30°, 5) anomalous refraction. 
The object is assumed to pass the meridian in the middle of a single observational run. We 
use 30 minutes of integration time or increase it such that the deconvolved maps look similar 
to the input model with chemical gradients. In contrast to modern interferometric facilities 
there will be no problem to achieve a good uv coverage with ALMA in a fraction of an hour. 

The simulated ALMA channel maps of HCO"^(4-3) at various spatial resolutions for the 
three adopted models and the inclination angle of 60° are presented in Fig. [3] (upper 3 rows). 
All significant features that are present in the "ideal" molecular emission spectra appear 
clearly at the resolution of 0.25" ("zoom-c" array configuration) and some are still visible 
at half the resolution ("zoom-d" configuration). The use of even lower resolution ("zoom- 
e") makes it difficult to disentangle the chemical and thermal features without extensive 
modeling. The DM Tau disk model with layered HCO^ abundances and 2D temperature 
gradient has the lowest line intensities and hence requires the longest observation of 2 hours, 
while for the models with uniform HCO^ abundances it can be as short as 30 minutes or 
less. 

In addition, we perform a similar analysis for a face-on orientation with i = 20°. At 
such low inclination channel maps for distinct disk models look similar and multi-line, multi- 
molecule analysis become a must. The corresponding V = 0.3 km s~^ channel of the HC0"'"(4- 
3) channel map for the disk model layered abundances is presented in Fig. [3] (bottom row). 
This channel map has a ring-like shape with two emission peaks located around the low- 
intensity midplane. The intensity in this channel is by a factor of 2 stronger than for the 
disk model with i = 60° and thus the observational time required to detect all these features 
is only 1 hour with the 0.25" beam size and less than 30 minutes at lower resolutions. All 
major features in the i = 20° channel map are resolved even with the compact "zoom-e" 
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array configuration. 



Given the importance of multi-line observations and ability of ALMA to simultaneously 
observe several transitions of several species in various frequency bands, we continue our 
study for other observationally important molecular tracers (CS, HON, CO isotopes) as well 
as a smaller disk with a radius of ~ 250 AU. All results are summarized in Table [H The 
channel maps for other transitions at millimeter wavelengths and for molecules located in the 
intermediate layer are similar to that of HCO"^, though their intensities are vastly different. 
We find that smaller disks are better studied at high frequencies of 400-700 GHz with 
moderately extended array configurations leading to smaller beam sizes than at millimeter 
wavelengths, albeit longer integration times are required for that. The contamination of the 
high-lying emission lines by the optically thick dust continuum from the disk inner regions 
and non-LTE excitation can be an issue for the analysis of these high-frequency data. 

The molecular lines also comprise kinematic information. High-quality ALMA spectra 
will allow to derive accurate dynamical ma sses of the central stars and thus calibrate pre- 



main-sequence evolutionary track models ( Simon et al.l 120001). d iscern gas motions along 
spiral arms in gravitationa lly unstable disks (IFromang et al.ll2004ar) . and trace flows of matter 



around a forming planet (jPaardekooper 200?! ). We do not pretend to analyze the ALMA 



requirements necessary for such studies as it would deserve a separate publication. 



5. Conclusions 

We showed that the ALMA interferometer will allow us to distinguish the effects of tem- 
perature gradients and chemical stratification in disks through molecular line observations, 
in particular for highly-inclined objects. We found that moderately extended array configu- 
rations (with baselines of < 1 km) and 0.5-10 hours of integration time will be necessary to 
pursue such observations. 

The authors are grateful to an anonymous referee for valuable comments and sugges- 
tions. This research has made use of NASA's Astrophysics Data System. 
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Fig. 1. — Distributions of the particle density (left panel), temperature (middle panel), and 
absolute abundances of HCO^ at 5 Myr (right panel) for the model most similar to the 
DM Tau disk. 
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Fig. 2. — (Left to right) The continuum-subtracted HCO"'"(4-3) synthetic map at the 
V = +0.77 kms~^ velocity channel for the three disk models: the model with chemical 
stratification and temperature gradients, the same model but with the uniform abundances, 
and the model with the uniform abundances and no vertical temperature gradient. The 
inclination angle is 60°. Intensity is given in units of radiative temperature (Kelvin). 
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Fig. 3. — (From first to third row) The same as in Fig. [2] but processed with the GILDAS 
ALMA simulator for the three array configurations and 64-antennas: "zoom-c" (~ 0.25" 
beam size, integration time is 2 hours), "zoom-b" (~ 0.5" beam size, integration time is 
0.5 hours), and "zoom-e" (~ 1" beam size, integration time is 0.5 hours). (Bottom row) The 
HCO"^(4-3) channel map at the V = +0.3 kms~^ velocity channel for the disk model with 
chemical gradients and the inclination angle of 20°. 



Table 1. Requirements for ALMA to study chemical and thermal gradients in disks 
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*Observing times are computed for the ALMA made of 64-antennas. With 50 antennas, these values should be 
increased by a factor of 1.5. 



